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ABSTRACT 



We have analyzed the properties of the Na D doublet lines at /1/15890,5896A in a large sample of 691 radio galaxies using the Sloan 
Digital Sky Survey (SDSS). These radio galaxies are resolved in the FIRST survey, have redshifts less that 0.2 and radio flux densities 
at 1.4 GHz higher than 40 mJy. The sample is complete within the main spectroscopic magnitude limits of the SDSS. Approximately 
1/3 of the sources show a significant excess (above that contributed by their stellar populations) of Na D absorption that can be 
robustly fitted with two Voigt profiles representing the Na D doublet. A further 1/6 of the sources show residual absorption, for which 
the fits were not well constrained though while ~50% of the sample show no significant residual absorption. The residual absorption 
is modestly blueshifted, typically by ~50 km s"', but the velocity dispersions are high, generally ~500 km s"'. Assuming that the 
size of the absorbing region is consistent with ~1 kpc for dust lanes in a sample of generally more powerful radio sources, assuming 
a continuous constant velocity flow (continuity equation), we estimate mass and energy outflow rates of about 10 Mq yr"' and few 
X 10*^ erg s"'. These rates are consistent with those in the literature based on Hi absorption line observations of radio galaxies. The 
energy required to power these outflows is on the order of 1-10% of the jet mechanical power and we conclude that the radio jet 
alone is sufficient. The mass and energy outflow rates are consistent with what is needed to heat/expel the mass returned by the stellar 
populations as well as the likely amount of gas from a cooling halo. This suggests that radio-loud AGN play a key role in energizing 
the outflow/heating phase of the feedback cycle. The deposition of the jet mechanical energy could be important for explaining the 
ensemble characteristics of massive early type galaxies in the local universe. 
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1. Introduction 
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The self-limiting cycle whereby the fueling of a supermassive 
black hole regulates both the rate at which it is fueled and the 
rate at which the surrounding spheroid grows - the feedback 
from active galactic nuclei (AGN feedback) - is engendering 
a substantial legacy in theoretical astrophysics. This cycle is 
invoked to explain many things, from why massive early type 
galaxies are "old, red, and dead" and thei r evolutionary behav- 
ior appears "anti-hierarc hical" (ICM; e.g. iThomas et al. Il2005t 
iBest et all l2005bl l2006l) . t he relationship betwe en the masses 
of black holes and bulges dTremaine et al. Il2002h . the "entropy 
floor" and lack of cool i ng flows in the inter- cluster medium (e.g. 
iFang & Hairnanll2008t iRaffertv et aDl2008h and pr eventing the 
coolin g of gas from the evolving stellar population dCiotti et al.l 
l2009l) . plus a seeming plethora of other characteristics of nearby 
galaxies, their halos, ICM and inter-galactic medium. 

Despite its potential importance, AGN feedback has re- 
mained in the realm of theoretical deus ex machina processes. 
This is mainly because it is difficult to obtain the necessary ob- 
servations that may tell us how AGN feedback regulates both 
the growth of the black hole and the galaxy. Here, we focus 
on radio sources because their mechanical energy output can 
be estimated, allowing us to gauge whether jets are a viable, 
though perhaps not t he only, m e chanism for creating the neces- 
sary feedback cvcle. lBest et al.l ( l2005bl) demonstrated that 30% 
of all nearby, early-type galaxies more luminous than L* in the 



SDSS have radio-loud AGN, and subsequently argued that this 
mec hanical energy o utput could balance cooling of the X-ray 
gas dBest et al.ll2006h . Powerful radio sources, based on ener- 
getic and (and in some cases) morphological arguments, ap- 
pear capable of driving energetic outflows of gas, which encom- 
pass a large number of ph ases including rnolecu lar ( Fisch er et al.l 
201 Oh. wa rm ionized (e.g.lFu & Stocktoiil2009l:iHolt et al.ll2008l : 
Nesvadba et al. 200M Fischer et alj|201 ll) . and warm neutral gas 
(e.g. lMorganti et a!]l2005allbll2007h . 



Absorption lines provide a robust way of probing one half of 
the feedback cycle. Because absorption lines are seen projected 
against the emission from the galaxy, negative radial velocities 
relative to systemic are an unambiguous sign of outflows. Most 
of the strong resonance lines of cosmically abundant ions are 
found in the UV (e.g., Morton 1991; Savage & Sembach 1996), 
and so must be studied from space, with, e.g., HST or FUSE. 
In the present study we have instead opted to exploit the vast 
data set of the SDSS to study a large sample of radio galaxies 
using the Na i doublet at /1/15890, 5896A, optical absorption lines 
that present another way of probing outflows from AGN. The 
ionization potential of Nai is only 5.1 eV, corresponding to a 
wavelength of A ~2420A, which is less than that of hydrogen. 
This implies that the photons that ionize Nai are in the near- 
UV and that sodium is mainly shielded from ionizing radiation 
by dust. These lines therefore primarily pr obe the dusty warm 
atomic phase and the cold molecular phase (ISpitzeiJll978l) . 
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For Na D lines to be observed, only relatively modest optical 
depths and H i column densities are required, which makes them 
a sensitive probe of the outflowing (or inflowing) neutral ISM. 
To observe the lines, the extinction must be suflicient for t>1 at 

2420A, which corresponds to Av=0.43 mag in the V-band for a 
ICard elli et al. ( 1989) extinction law and a selective extinction of 
Rv=3. 1 and to an H i column density of 8 x 10~ 



,20 ^^-2 



2. Sample and methodology 

The target sample is drawn from the seventh data release of 
the SDSS survey. All galaxies with spectroscopic redshifts were 
cross-matched with the NVSS and FIRST radio su rveys follow- 
ing the techniques outlined in iBest et al] (l2005ah . All sources 
have a 1 .4 GHz flux density (from FIRST and NVSS) exceeding 
40 mJy, i.e., more than an order of magnitu de less than those ob- 
served in Hi (iMor^anti et alJl2007l r2005allbl) . All galaxies have 
been visually inspected to ensure secure detections, and their ra- 
dio morphologies have been classified. The NaD sample consists 
of all 691 of these SDSS galaxies that have redshift z < 0.2 and 
an extended radio morphology at the k 5 arcsec resolution of the 
FIRST survey. We have focused on radio sources with extended 
morphologies beca use in these the jet mec hanical power can be 
robustly estimated JCavagnolo et al.lBOIOh . 

To estimate the contribution of any Nal in the interstellar 
medium of the host galaxies, we first had to remove the stel- 
lar contribution to the absorption lines by fitting each SDSS 
spectrum w ith a hnear combination of stellar population models 
taken from lBruzual & Charloll (l2003h using the p ublicly avail- 
able code STARLIGHT JCid Fernandes et al.l2005h . The regions 
of strong emission lines and those around the Na D lines were 
excluded from the fits. All spectra were well fitted and had re- 
duced ;i'^ values, suggesting the fits were highly significant. 

After removing the stellar continuum in the normalized spec- 
tra, we fitted both lines of the Na D doublet in a 100 A wide re- 
gion surrounding them, using a minimization routine and assum- 
ing the lines are repres ented by a Voig t profile with the appropri- 
ate atomic parameters (lMortonl2003l) . About half the sample did 
show significant Na D absorption line residuals, in fits that were 
robust based on their ;^'J values. However, because of the subtrac- 
tion of the stellar contribution of the lines, about 1/6 of the fits 
to the sources with significant residual absorption were not well 
constrained. This manifested itself most obviously when the core 
of the line was not deep compared to the overall rms noise level 
of the residual spectrum. We therefore adopted a threshold on the 
greatest depth of the lines to be used for further analysis of about 
5 times the rms of the residual spectrum. This is conservative, in 
that some of the sources with line cores weaker than this do show 
Na D absorption residuals. But given the reasonably large range 
of possible widths and offset velocities for such sources, we feel 
it is prudent to exclude them from the analysis (we refer to these 
sources as "unconstrained" for the sake of brevity). Furthermore, 
because the He i line at 5875. 6A overlaps with the blue wings of 
the (broad) Na D lines, our desire to obtain robust results made 
us exclude all 14 galaxies with strong He i emission. Thus, over- 
all, about 1/3 of all the sources had significant Na D residuals 
that had fits that were well constrained. The fit parameters are 
the line core depth, dispersion, and the velocity off'set relative 
to the systemic velocity (based on the redshifts determined from 
other stellar absorption lines in the spectrum). In all cases, the 
lines were sufficiently broad for the doublet not to be resolved 
into two distinct components. 
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Fig. 1. Distribution of fitted NaD absorption line offset velocities 
for our sample of radio sources. The classification of fit qualities 
is indicated in the figure legend. 
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Fig. 2. Distribution of NaD absorption line velocity dispersions 
for our sample of radio sources. The classification of fit qualities 
is indicated in the figure legend. 



Nai, our gas kinematics tracer, is a minor constituent of 
the warm neutral medium, most of which is in H i. To estimate 
the column density of H i relative to Na i, we require a Na de- 
pletion factor for atoms lost on grains, an ionization correc- 
tion and a Na/H abundance ratio. While we cannot determine 
these quantities directly, if we assume the clouds were similar to 
what is observed in the Milky Way, with a solar Na to H abun- 
dance ratio (-5.70 in the log), then the depletion on grains is 
likely to be abo ut a factor of a few and the ionization correc- 
tion about 100 (iPhillips et al.llI984l) . However, there is a weak 
trend of decreasing Na abundance with increasing column den- 
sity (Wakker & Mathis 2000) , and over the range of likely Hi 
column densities in our tar gets, the Na i column de nsities are a 
few hundred times lower (IWakker & Mathisl l2000l and refer- 
ence therein). Based on these estimates, we have adopted a ratio 
of ~460 to account for depletion and ionization correction. 

3. Results 

In our sample of nearly 700 radio-loud early-type galaxies at 
z<0.2, with 1.4 GHz flux densities above 40 mJy, about 1/3 
(260) of the sources show a significant detection of the Na D 
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doublet in absorption, whereas half (362) the sample do not. We 
probed column densities log Nnqi (cm"^)»11.5-13.5, which, not 
surprisingly, are about the level necessary to provide Tjust = 1 at 
2420 A (5.1 eV). Converting our estimate of Nai column den- 
sities to H I (using the Milky Way-like conversion factors given 
in § [2ll, log Nhi (cm~ ^ )~20-2 2 similar to what is probed in H i 
(iMor^anti et al.l l2003l l2005allbl) . The velocity offsets are rela- 
tively small, i -200 km s~', and even positive for a handful 
of sources (Fig.[TJ. The velocity dispersions, o-nhd, are typically 
-500-600 km s' (Fig.|2]i. 

Our goal is to understand the energy and mass outflow rates 
that may, or may not, be driven by the radio jet. Key to this analy- 
sis is the robustness of our fits. We therefore constructed a Monte 
Carlo simulation that creates a series of artificial lines with the 
same characteristics as the observed sample, which were then 
fitted in the same way as the data. For this, we simply assumed 
either (1) a constant value for both the covering fraction and 
the column density (both quantities were varied between sim- 
ulations) and varying signal-to-noise ratios, velocity offsets and 
dispersions according to a uniform distribution with ranges like 
those obtained from the data, or (2) that all quantities are uni- 
formly distributed, similar to the range of values observed. Our 
fits to the observed data show a trend between the covering frac- 
tion, C/, and the column density, Nnbd (Fig. ID, and we needed 
our robustness simulations to assess if this is an artifact of the 
fitting or if it has an astrophysical cause. What we found (espe- 
cially from the simulations where we kept the column density 
and covering fraction constant) was that we can reproduce the 
characteristics of the relationship between C/ and NNao- At high 
values for both Nnbd (well above 10'"^ cm~^) and covering frac- 
tions (above about a few tenths), the fits are robust, do not show 
a particular strong correlation between Cy and Nnbd, and repro- 
duce the input distribution of the dispersion and velocity offset. 
It is in the lower column and/or lower covering fraction regimes 
that we see this correlation, which is purely artificial. This is 
simply because the depth of the line is a function of C/, crNao 
and Nnbd, with different dependences if the lines are saturated 
or not. When the line is not (strongly) saturated, C/ and NNao 
are highly degenerate. 

The most robust quantities we found for a good recovery of 
our input distributions were the dispersion, velocity offsets, and 
the product C/ x Nnsd- The latter is important because the en- 
ergy and mass outflow estimates depend on the combination of 
these quantities, he input distribution of the dispersions is gen- 
erally recovered, and only the number of galaxies with high dis- 
persions are slightly underestimated (the recovered means were 
robust). The combination of C/ x NNao is recovered statistically 
to within a factor of ~2. For example, if we assume a constant 
combination of C/=0.5 and Nnbd = 3 x lO'-' cm"^, and uniform 
distributions for the dispersions, velocity offset and S/N, we find 
that the resulting product, C/NNao, has a mean of 1.5 x 10'^^ 
cm"- and that 68% of the data are between 0.5-3 x 10'^ cm"- 
- i.e., a majority of the values lie within a factor of about 2 (the 
results are shown in Fig.[3]i. 



4. Discussion and conclusions 

We can estimate the outflow rates of mass, M, and energy, E, in 
these galaxies by assuming that the continuity equation for mass 
and energy applies. For a spherical, continuous gas outflow with 
a covering fraction C/, Hi column density Nhi, wind opening 
angle Q, flowing through a radius r, and having a constant veloc- 




l,v„„ (cm--) 



Fig. 3. Relationship between the fitted column density of Na, 
NNaD, and the covering fraction, C/, of the best-fit model for 
the Na D lines. The color and size of each circle is related to 
the velocity offset (purple to red indicates -100 to -h200 km s"') 
and dispersion of the lines (lowest to highest indicates 200 to 
~1000 km s"') respectively. This trend is artificially induced be- 
cause the individual components of the doublet are not resolved 
and does not affect the overall results. The inset (magenta histro- 
gram) shows the distribution of Nnbd x C/ from the Monte- 
Carlo simulations, where the column density and covering frac- 
tions were held constant at 3x10'^ cm"^ and 0.5 respectively 
(see text for details). 



ity equal to that obse rved, the continuity equation implies (e.g. 
lHeckmanetai1l2000l) 

M~60 (r/kpc)(NH/3- 102icm-2)(5v/200 km/s)(Q/47r)Cf (1) 

£-8x10^' (r/kpc)(NH/3- 102icm-2)(5v/200 km/s)3(Q/47r)Cf (2) 

The observed frequency of significant Na D absorption is 
strongly influenced by the opening angle (because the back- 
ground against which the absorption line is seen is large), 
whereas C/ influences (mostly) the depth of the Une. In addi- 
tion, to obtain a more adequate representation of the terminal 
velocity of the Na D lines, we have assumed that <5v is not the 
measured offset velocity but rather the geometric combination 
of the offset velocity and the line dispersion. Given the range of 
values, t he terminal velocities are on the order of 1000 km s ' 
(see also lMorganti et al.ll2005ah . 

De-acceleration of the radio jet and shock heating of the 
inhomogeneous circum-nuclear gas transfers both momentum 
and energy from t he jet to the ambient ISM (e.g., the numeri- 
cal simulations of IWagner & Bicknellll201 11) . An energy-driven 
bubble develops, which rapidly e xpands and further accelerates 
additional gas and ISM clouds ([Sutherland & Bicknell 120071: 
IWagner & Bicknelll201 lb . Because we cannot directly measure 
the radii over which the energy-driven bubble is entraining and 
accelerating clouds in any individual source, we must estimate 
these radii in a more circuitous manner. 

Most of the mass that is accelerated is located where the gas 
and cloud density is the highest - within the inner ~kpc. Clouds 
in the ISM are accelerated very quickly (<10^ yr) to high ve- 
locities, and over relatively small radii only (<kpc), as shown 
in numerical simulations by, e.g., Wagner & Bicknell (201 1). Of 
course, while these clouds accelerate away from the AGN, they 
are heated, ionized and destroyed by mechanical and thermal in- 
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Fig. 4. Blue histograms show the energy {top) and mass {bottom) 
outflow rates estimated for our sources with robust NaD detec- 
tions, given a flow radius of 1 kpc for afl sources. The estimated 
jet mechanical power is represented by the aqua histogram in the 
top panel, which can be compare d with that of radio sou rces with 
Hi absorption detections (white: iMorganti et al.ll2005al) . The jet 
power estimates would be abou t an order of magnitude greater 
if we used the relationship from lCavagnolo et al.l(l2010l) . 



stabilities (iNakamura et al.ll2006l) . We have also argued that be- 
cause of the low ionization potential of Na, Na D absorption 
must be sampling dusty clouds in the warm atomic and cold 
molecular phases of the ISM. Therefore, we adopted 1 kpc for 
the radius the absorbing material is flowing through, this being 
the average angular s ize of the dust lanes in the 3CR radio sam- 
ple imaged by HST dde Kofi" et al.ll2000 h. Because the probable 
1 kpc size of the dust lanes is much smaller than the average 7 
kpc projected size of the SDSS spectroscopic fiber, this favors 
a small covering fraction for the Na D absorption and there- 
fore high column densities. The absorption is likely dominated 
by clouds that have high column densities (~10^'"^^ cm"^) but 
a low covering fraction (~10%) of the optical continuum light 
from the host galaxy. 

We find that the mass and energy outflow rates are about 10 
Mq yr"' and a few xlO''^ erg s"', respectively. Our sources have 
an average and median radio luminosity at 1.4 GHz of 10^^ W 
Hz"' and 10^^'^ W Hz"', respectively. The mechanical energy 
from the jet can be estimated using a scaling relation ship with 
the 1.4 GHz luminosity. Adopting the one given in iBest et al.l 
(l2006l and references therein) yields a typical mechanical en- 
ergy of our sources of ~2 x 10"*^ erg s"' (Fig.|4]i, while the re- 



lationship from lCavagnolo et al.l (1201 Oh gives estimates that are 
about one order of magnitude greater (this difference in energy 
estimates is mainly due to differing assumptions about the work 
done in inflating the X-ray cavities observed in clusters). Given 
the order-of-magnitude spirit of these estimates, it appears that 
about 1-10% of the jet mechanical energy is needed to power 
these outflows. This modest requirement implies that it is quite 
plausible for the jet itself to provide the required power A simi- 
lar conclusion based on X-ray cooling arguments was also drawn 
bv iBest et al.l (l2006l) . We do not find a correlation between the 
energy of the outflow and the jet mechanical energy. This is 
likely owing to the crudeness of our estimates but also to real 
astrophysical effects such as the coupling efficiency of the jet 
mechanical energy to the ISM, the relative direction of the jet 
compared to the distribution of the ISM, the range of masses of 
the warm neutral medium in these galaxies, etc. 

The mass outflow rates derived in this an alysis are simi- 
lar to those estimated for radio sources in Hi (iMor^anti et al.l 
12005 alibi) . This is not surprising given the similarity in approxi- 
mate jet mechanical energies of the two samples (Fig.|4]i. Over a 
plausible lifetime of 10^"*^ yr, an extended radio source will drive 
out/heat about lO*^"** M0 of gas and inject about 10^^"^^ erg of 
energy. This estimate can be compared with the rate at which 
the stellar populations return gas to the ISM of their galaxies 
and the total binding energy of the galaxy. A massive (stellar 
masses a: 10" Mp)) evolved (~10 Gyr) galax y will return about 
1 Mq yr"' (e.g. iMathews & Brighenti II2003I) . While we do not 
know how much gas the jet will remove from the potential of 
the galaxy, our estimated mass and energy outflow rates are sub- 
stantial compared to both the total amount of ISM in the typical 
elliptical, 10^ M0 (which is mostly in X-ray emitting gas), the 
rate of return from the stellar population, and amounts to virtu- 
ally the binding energy of a massive galaxy over the lifetime of 
the AGN. 

In massive galaxies jBest et al.l (l2005bl) found at the radio lu- 
minosities spanned by our sample that about a few to 10% have 
active radio sources. This either suggests that the lifetimes of 
radio-loud AGN are very long (which is unlikely) or that the 
duty cycle of AGN is high. If the lifetimes are as short as 10^ 
yr, this implies that the inactive phase lasts for only a few 100 
Myr, during which period it is likely that the host would accu- 
mulate an additional 10^"^ M0 of gas if most of it were caused 
by mass loss from stars within the galaxy. This amount is again 
similar to what t he radio source will heat during its lifetime (see 
iBest et al]l2006i for more details). From these rough estimates 
we conclude that the AGN is powering at least half of the nec- 
essary feedback cycle - the outflow/heating phase. More obser- 
vations and better modeling are needed to constrain the whole 
cycle of AGN feedback (e.g. Nesvadba et al. 2010). 
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